Resonant tunneling magneto resistance in coupled quantum wells 
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A three barrier resonant tunneling structure in which the two quantum wells are formed by a 
dilute magnetic semiconductor material (ZnMnSe) with a giant Zeeman splitting of the conduction 
band is theoretically investigated. Self-consistent numerical simulations of the structure predict 
giant magnetocurrent in the resonant bias regime as well as significant current spin polarization for 
a considerable range of applied biases. 



Semiconductor spintronics offers new functionalities to 
the existing electronics technology by combining charge 
and spin properties of the current carriers 0]. The goal 
of spintronic devices is to modulate charge current by 
changing the electron spin. In magnetic tunnel junctions, 
for example, magnetic electrodes are sandwiched around 
a nonmagnetic tunnel barrier, yielding large magnetocur- 
rents (relative current magnitudes for parallel and an- 
tiparallel orientations of the electrodes' magnetizations) 
In other spinteomcjieyices^ e^g., magnetic resonant 



tunneling diodes [2j, 11 m la M 121 l| , the transmission at 
a given voltage can depend strongly on the spin orienta- 
tion of the electrons at the Fermi level. Such diodes have 
been used as spin filters or spin detectors. 

Here we propose to use magnetic resonant tunneling 
diodes comprising coupled magnetic quantum wells and 
three nonmagnetic barriers, as spintronic devices offer- 
ing large magnetocurrents. Such devices combine the 
geometry of magnetic tunnel junctions with magnetic 
resonant tunneling diodes. Corresponding three-barrier 
nonm agne tic structures have been experimentally stud- 
ied 
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magnetic three barrier structure 
allows to establish parallel and antiparallel magnetiza- 
tion configurations and to observe magnetocurrent. The 
origin of the magnetization here is not essential. The 
quantum wells can be either dilute magnetic semicon- 
ductors 0|, which have giant g- factors at temperatures 
up to 100 K, or they could be ferromagnetic quantum 
wells 0, 0, 0] in which the splitting is given by the 
exchange field. Performing realistic sclfconsistent calcu- 
lations of the I-V characteristics we indeed predict large 
magnetocurrents (say, 10000% for a moderate spin split- 
ting of 10 meV) at resonant voltages. These large values 
appear because of the suppression of the resonant tun- 
neling due to the energy mismatch for the resonant levels 
of the coupled wells in the antiparallel configuration. Al- 
though the predicted effect is rather robust, we provide 
hints on how to further tune the heterostructure param- 
eters to maximize the magnetocurrent. 

We investigate a three barrier semiconductor het- 
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FIG. 1: (Color online) Schematic conduction band profile of 
an all II-IV semiconductor three barrier heterostructure. The 
quasibound states in the quantum wells are Zeeman-split by 
an external magnetic field. 



erostructure based on a Zn(Be,Mn)Se material system. 
The conduction band profile at zero bias is schematically 
shown in Fig. ^ The barriers are formed by Be-doped 
ZnSe, whereas the quantum wells (QWs) are made of 
the dilute magnetic semiconductor (DMS) material Zn- 
MnSe. Such a system is an extension of a two barrier 
structure already experimentally inve stig ated in Ref. 0] . 
Following a previous theoretical work [ljj on two-barriers 
structures, we assume the barrier height to be 92 meV, 
which is about 23% of the band gap difference between 
ZnSe and ZnBeSe. In DMSs an external magnetic field 
B causes a giant Zeeman energy splitting AE into spin 
up and spin down electron states, which can be expressed 
by the modified Brillouin function B s |8(: 



AE = x o sN asB s (g M nSIJ-BB/kBT eS ). 



(1) 



Here, x D s is the effective concentration of Mn ions, 
Noa = 0.26 eV is the sp—d exchange constant for conduc- 
tion electrons, s = 5/2 is the Mn spin, gun — 2-00 is the 
Mn g- factor, {1b labels the Bohr magneton, ks denotes 
Boltzmanns constant, and T e g is an effective tempera- 
ture. At low temperatures (all simulations are performed 
at a lattice temperature of T= 4.2 K), Mn concetrations 
of about 8%, and practical magnetic fields of 1-2 T, the 
giant electron g-factor in the DMS leads to an energy 
splitting of the conduction band of about 10 meV |2|. 
The whole structure is considered to be sandwiched be- 
tween two leads consisting of n-doped ZeSe (n = 10 18 



2 



cm~ 3 ) of 10 nm width. Similar to experiments 0], we in- 
clude 5 nm thick undoped ZnSe buffer layers between the 
leads and the active structure, which causes an upward 
band bending at zero bias. 

Following the classic treatments of (two barriers) res- 
onant tunneling diodes 0, 0, |2jj , we assume coherent 
transport throughout the whole active region. This is a 
reasonable assumption also for our three barriers mag- 
netic diode, since signatures of coherence, namely the 
splitting of the well states into bonding and antibond- 
ing states, have been observed experimentally in such 
structures Q. We calculate the spin dependent current 
flow by numerically solving the conduction band effective 
mass Schrodinger equation taking into account the spin 
dependent potential energy, 

Ui{z) -e<j)(z) + aAE{z). (2) 

Here, Ui is the intrinsic conduction band profile, e is 
the elementary charge and a = ±1/2, (T, I) labels the 
spin quantum number. The electrostatic potential <fi is 
obtained from the Poisson equation, which is solved to- 
gether with the Schrodinger equation in a self-consistent 
way. The current density of electrons with spin a is cal- 
culated by 

ja= (2^)3 / dSkv M E hEt)\f(E)-f(E + eV a )], (3) 

where Ei and E t are, respectively, the longitudinal and 
transverse component of the electron total energy E, 
V a denotes the applied voltage, T a (Ei,E t ) is the elec- 
tron transmission function, v z labels the longitudinal 
component of the electron group velocity, and f(E) = 
1/[1 + exp(_E - Ef)/k B T] is the Fermi function with 
the fermi energy Ef. Assuming parabolic bands and 
using the same effective mass m,/mo = 0.145 (with 
mo denoting the free electron mass) for all layers of 
the heterostructure, the transmission function only de- 
pends on the longitudinal energy, T{Ei 1 E t ) — T{E{). 
This allows to reduce Eq. to the Tsu-Esaki formula 
|2l| . The current spin polarization is then determined by 

We assume that the magnetization of the first QW 
is fixed, whereas the second is 'soft', which means that 
it is sensitive to local changes of an external magnetic 
field. Hence, a magnetocurrent MC can be defined as 
the relative difference of the current I for parallel (P) 
and antiparallel (AP) alignment of the magnetization of 
the two quantum wells, MC = (Ip — Iap) I Iap ■ The pro- 
posed structure aims in producing very high MCs based 
on the following idea of operation. High resonant tunnel- 
ing throughout the whole structure is possible if the mid- 
dle barrier is thin enough and if two quasi-bound states of 
same spin of the adjacent QWs are aligned energetically. 
Such resonant condition for the lowest energy states in 
the case of P magnetization is illustrated in the inset of 
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FIG. 2: (Color online) The spin resolved energy-dependent 
transmission function at T = 4.2 K for parallel (P) and an- 
tiparallel (AP) magnetization in the case of a given Zeeman 
splitting of AE = 10 meV and a resonant applied voltage 
of V a = 13 meV. The inset shows a contour plot of the local 
density of states versus energy and growth direction z at reso- 
nance (V a = 13 mV) for parallel magnetization and AE = 10 
meV. The solid line indicates the self-consistent conduction 
band profile. The higher quasibound states (at about 80 and 
90 meV) are not in resonance. 



Fig. |21 The inset shows the local density of states of the 
conduction electrons and clearly demonstrates the Zee- 
man splitting of the quasibound states into spin up and 
spin down states. Here, the orientation of the magnetic 
field is assumed such that the spin down (up) energy lev- 
els are shifted downwards (upwards) in energy. Due to 
the interaction of the two QWs the resonant energy lev- 
els are further split into bonding and antibonding states. 
The splitting is observable in the local density of states 
if it is greater than the natural energy broadening of the 
quasibound state. For our structure the middle barrier is 
too thick to resolve this additional splitting in the plot of 
the local density of states. The inset of Fig. [21 also shows 
higher quasibound states at about 80 and 90 meV, which 
are however not in resonance. 

In the case of P magnetization the spin up and down 
quasibound states are equally Zeeman shifted in both 
QWs, whereas for AP alignment they are shifted energet- 
ically in opposite directions. Assuming QWs of the same 
width, i.e, with the same quasi-bound energy spectrum, 
there are hence at equilibrium open resonant conduction 
channels for the P alignment, whereas the transmission 
is blocked for the AP configuration. However, when a 
small finite voltage is applied to the structure the res- 
onant channels of the P magnetization are 'destroyed', 
since the energy levels of the second QW are shifted more 
deeply to lower energies than those of the first QW by the 
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FIG. 3: (Color online) Current- Voltage characteristics of the 
structure for parallel (P) and antiparallel (AP) magnetization 
at the temperature T = 4.2 K and the Zeeman splitting of 
AE = 10 meV. The inset displays the current polarization as 
a function of the applied voltage. 

applied bias. To overcome this shortcoming, we propose 
to use asymmetric QWs. Here, we take the second QW to 
be thinner than the first one. This gives rise to a higher 
ground state energy in the second QW and the resonant 
condition is therefore adjusted at a finite voltage, leading 
to high currents. 

In order to maximize the MC we pursue the following 
strategy for determining the different layer widths. At 
low temperatures the current is given, in good approxi- 
mation, by the quadrature of the transmission function 
over the energy window A/i = [Ef ~ eV a , Ef]. To obtain 
high currents for the P alignment, the layer widths should 
be chosen such, that the resonant tunneling condition, 
T(E) s=s 1, is fulfilled at a finite voltage for energies be- 
longing to the interval A/x. On the other hand, for the AP 
configuration the transmission function should be made 
as small as possible in the energy window Afj,. To meet 
both demands at the same time, we use a relatively thin 
middle barrier, which effectively controls the coupling of 
the QWs. Figure |21 shows the spin-resolved transmission 
function versus energy for the resonant voltage V a = 13 
mV. The double peak structure of the transmission func- 
tion for the AP configuration corresponds to the low- 
est quasibound state in both QWs, which have in that 
case different energies leading to two 'half-resonances'. 
The transmission for the AP alignment can be strongly 
hampered by choosing thick QWs. The variation of the 
thickness of the first and third barrier barely influences 
the MC, since the transmission is then either reduced or 
increased for both P and AP magnetization. By chang- 
ing the buffer layer thicknesses one can tune the amount 
of band bending and hence, the relative position of the 
quasi-bound energies to the Fermi level. The dimensions 
of the structure finally used in our numerical calculations 
are indicated in Fig. ^ 

The obtained current- voltage characteristics for P and 



FIG. 4: (Color online) Magnetocurrent (MC) as a function of 
applied voltage for different Zeeman splittings AE. 



AP magnetization and a fixed Zeeman splitting of AE — 
10 meV are displayed in Fig. [3] In the case of P magne- 
tization, the lowest quasibound energy levels of the adja- 
cent QWs are already aligned at a small voltage of about 
13 mV, whereas for the AP configuration a much higher 
voltage is necessary to obtain resonant tunneling. The 
current spin polarization is plotted in the inset of Fig|3] 
Since the spin down state has a lower energy than the 
spin up state, the current for P magnetization is almost 
all spin down polarized at low voltages. For higher volt- 
ages also spin up states contribute to the current thereby 
diminishing the polarization. In the case of AP magne- 
tization, spin up current can flow at low voltages due to 
the lowest spin up state in the second QW. By inreas- 
ing the voltage this state gets off resonance, which leads 
again to significant spin down polarization. The obtained 
MCs for different Zeeman energy splittings are shown in 
FigEl The transmission for AP magnetization is strongly 
reduced by increasing the energy splitting AE. Hence, 
the current for AP magnetization becomes very small at 
the peak voltage of the P configuration and our simula- 
tions reveal very high MCs up to 800 for reasonable spin 
splitting. Simulations performed at a higher temperature 
T = 100 K show that the MC is reduced to about 25% of 
its value at T — 4.2 K. Although the temperature effect 
is quite large, the MC remains significant. 

To summarize, we have numerically investigated an 
all semiconductor three barrier resonant tunneling struc- 
ture, comprising two quantum wells made of a magnetic 
material. Our simulations predict very high magnetocur- 
rents demonstrating the device potential of the proposed 
structure. Further investigations of an in-plane (perpen- 
dicular to the growth direction of the heterostructure) 
resistance of the two coupled magnetic QWs promise to 
reveal interesting effects, since resistance resonances for 
such a current-in-plane configuration have already been 
demonstrated for nonmagnetic structures |22| . 
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